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f^ ' Recent research on the buildup of rocks from smaU dust grains has reaffirmed 

0\ I that grain growth in protoplanetary disks should occur quickly. Calculation 

Qs , of growth rates have been made for a variety of growth processes and gener- 

■"■^ — , ■ ally predict high probabilities of sticking in low-velocity collisions that may 

i-Q ' be brought about in a number of ways in protoplanetary disks. Laboratory 

C*h' experiments have measured sticking coefficients for some materials largely 

confirming the calculations. Although the detailed velocity fields of disks 
;**! . are not well understood, many of the important processes leading to particle 

collisions and grain growth have been studied theoretically and demonstrate 
likely paths by which dust is assembled into planets. Calculations of the ra- 
diative properties of particles with various size distributions show that large 
particles should produce observable changes in the spectral energy distribu- 
tions of disks. Changes of the sort predicted are, in fact, observed, but their 
interpretation is ambiguous; there are other ways to produce the observed 
JH , changes that do not require grain growth, so the evidence is currently incon- 

elusive. The major uncertainties can be overcome with the next generation 
of millimeter-wave interferometers, and it seems likely that a firm case for 
grain growth could be established within a decade. 



I. Introduction 

Planets are built from tiny interstellar dust grains by the steady ac- 
cumulation of the smaller particles, at first through sticking but later 
aided by gravity between larger planetesimals, and subsequently by the 
accretion of gas onto the rocky cores. Such growth should occur natu- 
rally in the circumstellar disks created when young stars are born. The 
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disks provide protected " wombs" for the development of the large bod- 
ies during the course of a few million years or more. Although disks are 
not the only places where planets might have grown, they are theoret- 
ically well suited for planet formation and observed to be ubiquitous, 
so that they are the clear favorite for planetary birthsites at the time 
of this review. A demonstration of the growth of dust into large bodies 
is the one open issue that remains to demonstrate that circumstellar 
disks - hundreds of examples of which are routinely studied - do spawn 
planetary systems. Such a demonstration would also imply that plan- 
etary systems are common features of nearby stars, just as disks are 
commonly found around young stellar objects. 

In the earliest stages, particles are built up by coagulation and 
sticking as they bump into one another in the disks. In about 10* yr, 
particles can grow from micron-size to meter-size through a variety of 
coUisional processes: Brownian motion, settling, turbulence, and ra- 
dial migration, the dominant mode depending on physical conditions 
within the disk. After another 10"* years or so, the small bodies can 
attain sizes of order kilometers. Beyond 1 km, gravitational attrac- 
tion causes these planetesimals to grow by pairwise collisions and later 
by runaway growth until some planets will be present after 1 Myr 
(WeidenschiUing 1988; Weidenschilling & Cuzzi 1993; Lissauer 1993). 
Accretion of gas should create gaps in the disks that may present a nat- 
ural limit to the growth of planets over very long times. Figure 1 shows 
the three stages schematically. The timescales are uncertain but still 
well enough constrained to be shorter than the observed lifetimes of 
disks around young stars (e.g. Strom et al. 1993). If the theory is even 
approximately correct, there should be large bodies within currently 
observed circumstellar disks. 

The buildup of a small, rocky planet is longer than the lifetime 
of any observer by several orders of magnitude, so it is not possible 
to watch this evolution. Instead, we must be content to compare the 
particle properties of disks of different ages and infer the growth of large 
bodies from trends in the size distribution with age. Even the indirect 
detection of planetary-mass bodies in circumstellar disks is difficult 
when the disk mass is greater than the planetary mass, but there is 
hope that the indirect effects of larger bodies will be observable in a 
decade or sooner. 

The theory of grain coagulation delineates the time scales over 
which to study disks and provides some limits on what is plausibly 
observed. The theory is checked by laboratory experiments for small 
scale growth, where major bottlenecks might occur. It is possible to 
create tiny grains similar to interstellar dust in the laboratory under 
physical regimes like those thought to occur in disks. The growth rates 
of these particles can be checked against the theory to ensure that the 
physics of grain growth is well understood. Experimentally measured 
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growth rates for small particles largely verify the calculations and lead 
us to believe that the growth of the smallest particles is reasonably well 
understood. It is difficult to verify the growth of larger bodies. It will, 
therefore, be important to detect some large bodies near other stars to 
check the calculations of early planet evolution. 
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Figure 1. This figure depicts the sequence of events leading to planet for- 
mation: grain growth through particle-particle interactions, growth of plan- 
etesimals by two-body interactions aided by gravity, and the accretion of gas 
via gravitational attraction. 

There are enough known disks around nearby young stars to pro- 
vide good samples for detection of grain growth. There are now meth- 
ods for discerning the growth of dust grains to pebbles and rocks 
(Miyake & Nakagawa 1993) via the submillimeter opacities. These 
methods have had limited application to date but give interesting hints 
of grain growth. 

This review discusses several aspects of this problem, including the 
physics of the growth process, the observable effects on radiation from 
collections of the particles, and the extant observations designed to see 
evidence of grain growth. There has been progress on all these fronts 
since the last Protostars and Planets meeting in 1990, spurred by an 
increasing interest in extra-solar planet formation that is driven by 
detections of extra-solar planets. Most of this progress has been in lab- 
oratory work and new calculations of the properties of small particles 
as they pertain to coagulation and growth. There were already tan- 
talizing hints at the time of the last conference, Protostars & Planets 
III that particle growth had been detected in circumstellar disks; a few 
new observations attempt to clarify these hints. We will concentrate 
on progress made since PPIII. 
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II. Grain evolution in disks 

A. Structure of protoplanetary disks 

The idealized infinitely thin, continuous disk assumed in early theories 
(e.g. Lynden-BeU & Pringle 1974; Pringle 1981; Adams & Shu 1986) 
has given way to a disk with vertical and radial structure and with 
gaps or holes needed to account for various observations (Bodenheimer 
1995; Papaloizou & Lin 1995). Stellar radiation provides at least min- 
imal heating that is often enhanced by accretion energy as material 
from the disk falls onto the star. Disks heated mainly by accretion 
are referred to as active, whereas those heated entirely by the stars 
are called passive. The spectral energy distributions of most simple 
disks can be fitted reasonably well assuming they are flared (Kenyon & 
Hartmann 1987), centrally heated from radiation near the star (Chiang 
& Goldrcich 1997) and have radially decreasing temperatures (Adams, 
Lada, and Shu 1988). There arc striking exceptions to this picture - 
HL Tau, for example, and almost all of the "flat-spectrum" sources; 
these presumably have a number of different components - dense, disk- 
like regions, outflow cavities, the surrounding molecular cloud cores - 
that contribute to the SEDs and complicate the analysis (Calvet et al. 
1994) 

Disks usually have outer radii between a few and a few hundred 
AU (Appenzeller et al. 1984; Edwards et al. 1987; Lay ct al. 1997) 
often with sharp outer boundaries (O'Dell & Wenn 1994; McCaughrean 
& O'Dell 1996). Young disks are thought to be optically thick at all 
wavelengths shortward of 100 /im. Typical disk masses are between 
0.001 and 0.1 M0, although there is some uncertainty about this range 
(Adams etal. 1990; Beckwithet al. 1990; Henning et al. 1993; Andre & 
Montmcrle 1994; Osterloh & Beckwith 1996). These properties accord 
well with those assumed for the disk in the early solar system from 
which our planets were born (Safronov 1969; Hayashi 1981; Cameron 
1988). 

Disk structure is important to the study of particle growth in two 
ways. First, the presence of convection, turbulence eddies, and other 
velocity gradients can strongly affect the growth of particles, since the 
growth rates depend on the relative particle velocities. Second, obser- 
vations of disks must either take disk structure into account for their 
interpretation or else be designed to be relatively unaffected by the 
structure. These effects are discussed in the sections to follow. There 
is enough freedom in models of disks to produce many kinds of local 
structure, and important effects such as turbulence eddies, sedimenta- 
tion, and radial migration of matter are currently beyond observation. 
They are also a matter of debate among theorists. While several insta- 
bility mechanisms might in principle maintain nebular turbulence, the 
presence of turbulence and the dominance of any energy source has yet 
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to be demonstrated definitively (Cabot et al. 1987; Balbus & Hawley 
1991a, b; Dubrulle 1992, 1993; Tscharnuter & Boss 1993; Nakamoto & 
Nakagawa 1994, 1995; Kley et al. 1993; Stone & Balbus 1996; Klahr et 
al. 1998). As a result, there is considerable uncertainty about the ac- 
tual evolution of particles in disk theories, and a wide range of physical 
conditions needs to be considered in theories and laboratory studies of 
grain growth relevant to planet formation. 

B. Physical processes leading to grain growth 

In the outer regions of an accretion disk, most of the interstellar dust 
remains unaltered during passage through the weak shock front when 
the disk is first created. Only volatile ice mantles around the refractory 
grain cores partly sublimate, but water ice efficiently forms when the 
water vapor cools behind the shock (Lunine et al. 1991, Neufeld and 
Hollenbach 1994). Pollack et al. (1994) summarize the expected grain 
components, their optical properties, and sublimation temperatures im- 
mediately after disk formation. The main species are: (1) olivines 
(Mga;Fe2_j:Si04 with x==1.4), (2) orthopyroxenes (Mga;Fei_a;Si03 with 
x=0.7), (3) quartz (Si02), (4) metallic iron (Fe), (5) troilite (FeS), (6) 
refractory and volatile organics, and (7) water ice. Henning and Stog- 
nienko (1996) update these calculations by taking into account more 
recent determinations of optical constants and a possible fluffy struc- 
ture of the grains. They also found that the grain opacities depend on 
the fraction of iron that is in iron grains or silicates. 

If the grains are transported inwards during the viscous evolution, 
the higher temperatures will anneal and eventually destroy them (see, 
e.g., Finocchi et al. 1997, Gail 1998). Gail (1998) assumes chemi- 
cal equilibrium and concludes that crystalline silicates, iron, and alu- 
minium compounds are the most important grain components in the 
warmer parts of the disk. Such equilibrium condensation calculations 
have been extensively performed (see, e.g., Grossman 1972, Cameron 
and Fegley 1982, Yoneda and Grossman 1995). The results are uncer- 
tain, however, because a good understanding of the chemical processes 
during grain condensation, growth, and destruction requires a kinetic 
approach which cannot yet be done owing to unknown reaction rates 
for oxygen-rich chemistry (see also Prinn 1993). 

The condensation of ice onto refractory cores in the outer regions 
of the disk is important for the growth of submicron-sized grains but is 
irrelevant to growth from micron-sized grains to meter- and kilometer- 
sized planetesimals (e.g. Preibisch et al. 1993). But the surface struc- 
ture of the dust grains may be considerably modified by repeated con- 
densation and sublimation caused by transport processes or heating 
events. Such "fiares" were also discussed as a possible mechanism for 
the production of chondrules found in meteorites (see Hewins et al. 
1996). 
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For much of the last 30 years, planetesimal formation was thought 
to be the resuh of a gravitational instability in a dense dust layer lead- 
ing to fragmentation into kilometer-sized bodies (Safronov 1969, Gol- 
dreich and Ward 1973). Weidenschilling (1980) was the first who noted 
the difficulty with the particle layer instability. Cuzzi et al. (1993) and 
Weidenschilling (1995) demonstrated that such an instability is unlikely 
to occur. Particles must grow to at least meter size before they can 
decouple from shear-induced turbulence and slow to the low velocities 
required for the instability to occur. The radial velocity dispersion in- 
duced by drag will further delay the onset of instability until the mean 
size is in the range of 10-100 m. Thus, the formation of planetesimals 
has to be explained by the collisional aggregation of particles. 

It is difhcult to calculate this process from first principles, be- 
cause there are complex feedback mechanisms in the two-fluid system 
whose relative importance is not well known: aerodynamical drag, ther- 
mal structure coupled to changing dust opacities, charge state of the 
medium and coupling to magnetic fields, for example. It is imperative 
to understand particle growth from collisions between particles, how- 
ever, if we want to estimate the likelihood of creating planetesimals and 
planets. 

Particles collide because they are moving relative to one another. 
The sources of relative velocities in laminar nebulae are Brownian mo- 
tion, vth, sedimentation, «„, and radial drift, Vr- In turbulent systems, 
turbulence-driven motion, vtur, can additionally be of importance for 
grains/bodies of certain sizes. The equations for the relative velocities 
are: 



Awt/i = y/SkT/niJ, with ^ = {mim2) / {mi + 1112) , (1) 

AVs=n'^Z \Tfl-Tf2\, (2) 

Aw^ = 2f7At/ \Tfi-Tf2\, (3) 

Avtur -Ujts IT/1-T/2I for Tfi,Tf2<ts, (4) 



< Avtur >= i-SSui ^Tf/ti for ti > Tfi,Tf2 > ts, (5) 

(WeidenschiUing 1977, Volk et al. 1980, Mizuno et al. 1988). Here, 
51 is the Keplerian angular frequency; T the temperature; mi and ?7i2 
the particle masses; r/i and r/2 the friction times (see IIIA for a more 
detailed discussion), r/ = max(T/i, r/2); z the vertical distance from 
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the midplane; Ay the velocity difference between Keplerian velocity 
and rotational gas velocity; Us and tg the velocity and turn-over time of 
the smallest turbulent eddy; ui and ti the velocity and turn-over time 
of the largest turbulent eddy. In the case of turbulence-driven motion, 
we have distinguished between two cases of the coupling of the grains 
to the gas (< Avtur > is the rms value of the relative velocities). 

A comparison of equations (1) through (5) shows that the par- 
ticles with equal friction times must couple to the turbulence (second 
case, eq. 5) to obtain substantial relative velocities; otherwise, only the 
Brownian motion is important. The collisions caused by Brownian mo- 
tion are relevant only for small (about micron-sized) grains because of 
the inverse mass dependence. On the other hand, an efficient coupling 
to the turbulence in a disk with q:=0.01, say, is possible for mm-sized 
grains to m-sized bodies, calculated for the disk model of Bell et al. 
(1997) at 10 AU radius and considering the friction time for compact 
particles. Grains with very small friction times essentially comovc with 
the gas. One needs friction times larger than the turn-over time of the 
smallest eddies for the start of an efficient growth by turbulence. In a 
laminar nebula, a dispersion in friction times is always required for an 
efficient growth. 

C. Homogeneous growth regime 

Most models for coUisional grain growth in protoplanetary accretion 
disks assume a locally homogeneous distribution of particles. There 
are classically two types of models: radial one-dimensional (ID) models 
that average the density distribution vertically through the disk (sec, 
e.g., Schmitt et al. 1997), and vertical ID models for the coagulation 
of particles by settling to the midplane at a given radius (see, e.g., 
Weidenschilling and Cuzzi 1993). 

Several authors have studied the ID radial case each incorporating 
more detail into the calculations (see, e.g., MorfiU 1988; Mizuno et al. 
1988; Mizuno 1989). Ruden and Pollack (1991) treated the dynamical 
evolution of a protoplanetary disk under the hypothesis that the only 
source for the turbulent viscosity is thermal convection. They discussed 
grain settling and coagulation in the outer solar system when turbu- 
lence ceases. Schmitt et al. (1997) coupled the disk and dust evolution. 
They distinguished three phases during the first 100 years of disk evolu- 
tion: an early phase in which the smallest particles disappeared quickly 
due to Brownian-motion driven coagulation and produced a relatively 
narrow mass distribution; an intermediate phase of self-similar growth 
due to turbulence where the size distribution can be described by scal- 
ing laws; and a late phase where the most massive particles decouple 
from the gas and drift motions become an important source of relative 
velocities. Stepinski and Valageas (1996) consider the global evolution 
of single-sized non-coagulating particles in radial models with a be- 
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tween 0.01 and 0.001. They conclude that particles larger than 0.1 cm 
are not really entrained in the gas. They also studied the formation 
and radial distribution of km-sized icy planetesimals assuming that the 
size distribution is always narrowly peaked around a mean value for a 
given radial location and time (Stepinski and Valageas 1997). 

A shortcoming of radial models is the inability to treat the forma- 
tion of a dust subdisk due to settling. Weidenschilling (1980), Nak- 
agawa et al. (1981, 1986), and Weidenschilling and Cuzzi (1993) de- 
veloped numerical ID models for settling with coagulation. The main 
steps of the evolution are: (1) Micron-sized grains uniformly suspended 
in the gas grow by Brownian motion, (2) grain aggregates begin to set- 
tle and sweep up smaller ones in a kind of runaway growth, and (3) the 
aggregates grow to centimeter-sizes and create a dense dust layer (Wei- 
denschilling 1997). Weidenschilling (1997) treatede the growth from 
microscopic grains to cometesimals at 30 AU in detail. The depen- 
dence of drag-induced velocities on particle size produced a stage in 
which most of the mass is concentrated in a narrow size range between 
10 and 100 m. DubruUe et al. (1995) found that the turbulence does 
not influence the timescale of sedimentation, but merely determines 
the equilibrium scale height of the dust subdisk. The smallest grains 
that can experience some kind of sedimentation in the presence of tur- 
bulence have sizes of the order of 150a cm, resulting in a size limit of 
about 1cm for q;=0.01. Sekiya (1998) showed that for mm-sized ag- 
gregates even very weak turbulence stirs up the dust aggregates. This 
calculation considered only locally produced turbulence corresponding 
to the calculations by Cuzzi et al. (1993) and Champney et al. (1995), 
not an a-disk. 

Radial transport becomes important over the timescale of settling. 
The results of the ID investigations show that consistent 2D/3D models 
or at least 1-l-lD models are necessary for an adequate description of 
the dust dynamics. 

D. Inhomogeneous growth regime 

In the classical models, turbulence is treated only statistically as a 
source of random velocities important for diffusion and the calcula- 
tion of turbulence- induced relative velocities. These models assume 
that particles remain uniformly distributed in a turbulent velocity field. 
This prediction is only valid if the typical length scales are larger than 
the eddy scales. 

Squires and Eaton (1991) were the first to show that turbulence 
can lead to particle concentration. Cuzzi et al. (1996) calculated this 
particle concentration for a protoplanetary disk with weak Kolmogorov- 
type tubulence. Concentrations of 10^ to lO'' occur for particles with 
radii between 0.1 and 1 mm in convergence zones between eddies. The 
particle trapping mechanism may have dramatic consequences for the 
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formation of chondnilcs due to its size-selective nature (see, Cuzzi et 
al. 1996). 

Klahr and Hcnning (1997) showed that in certain vortical flows 
particles are trapped and concentrated in the gas flow. Slowly rotating 
eddies and a gradient of the vertical component of the gravitational 
force produce this trapping. The gradient causes particles in the upper 
half of an eddy to fall faster than in the lower half. If this concentration 
is stronger than the dispersion of grains due to centrifugal forces, the 
particles are concentrated in the interior of the eddies. The strongest 
effect was observed for millimeter-sized particles; they can be concen- 
trated by a factor of 100 within 100 years. Three-dimensional calcula- 
tions by Klahr et al. (1998) indicate that the particle concentration 
mechanism is not just the result of 2D approximations. Figure 2 shows 
the results of a calculation in which an initially homogeneous distribu- 
tion becomes inhomogencous after 160 years. However, we should note 
that the possibility to concentrate or to disperse small dust particles 
depends on the structure of the flow pattern (3D disk turbulence vs. 
slow and long-lived vortical circulation) and the stability of the vortices 
(Adams and Watkins 1995; Klahr and Henning 1997; Klahr et al. 1998; 
Brandenburg and Hodgson 1998). 
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Figure 2. The left hand portion shows an initially homogeneous distribution 
of grains that becomes quite inhomogeneous after 160 years in the right hand 
portion because of trapping in eddies (Klahr et al. 1998) . 



Tanga et al. (1996; see also the earlier related study by Barge and 
Sommeria 1995) considered one more possibility of particle concentra- 
tion via vorticity. They introduced an ad hoc system of regular vortices 
driven by differential rotation of the disk with rotational axes parallel 
to its angular momentum. In this case, the concentration mechanism 
is based on Coriolis forces and much larger particles would be concen- 
trated compared with the cases discussed before. 
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III. How well is particle growth understood? 

A. Aerodynamical behaviour of particles in a dilute gas 

The interaction of dust particles and small bodies with sizes below 
1 km is determined by their aerodynamical coupling to the dilute gas 
of the protoplanetary accretion disk and not by gravitational forces. 
This coupling can be characterized by a friction (or response) time 
Tf which for dust particles and aggregates is much smaller than the 
orbital period; i.e. the particles are strongly coupled to the gas. For 
the conditions in protoplanetary disks, the coupling of submicron- to 
centimeter-sized particles to the gas is dominated by collisions of grains 
with individual molecules. In this regime, the molecular mean free path, 
Imfp, is larger than the characteristic particle size, a. The friction time 
of a spherical dust particle of mass, m, and geometrical cross section, 
tTg, is given by: 

Here SP is the momentum transfer efficiency factor and the quantities 
Pg and Vg denote the gas density and the mean thermal velocity of 
the molecules, respectively. The material parameter 6P is always close 
to 1. Experiments by Hutchins et al. (1995) and Blum ct al. (1996) 
verified the validity of (6) and showed that a single value 6P— 1.11 ± 
0.17 describes both spherical Si02 grains and small aggregates. These 
experiments confirmed the results of numerical Monte-Carlo simula- 
tions on the aerodynamic properties of fractal aggregates (Mcakin et 
al. 1989). 

The strength of gas-grain coupling depends on the ratio of mass 
to cross section, so t/ depends on whether the particles are compact 
- spheroids, say - or fluffy, open structures (Ossenkopf 1993). If the 
aggregates are fractals, we can introduce two fractal dimensions D„i 
and Da^: 

m ex i?^" and (jg ex R^" (7) 

and show that r/ oc '^(^^n.-D^) ^ ^n 3-dimcnsional systems, 1 < D^ < 
Dm < 3. Note that Dm = 3 and D^ = 2 for compact particles, so that 

Tf OC R. 

B. Aggregation process 

Ballistic cluster-cluster aggregation, denoted BCCA, is the process in 
which clusters that collide stick together without restructuring and 
with sticking probability 1. This process leads to the production of 
open, fluffy aggregates. Numerical simulations of BCCA yield Dm 
values between 1.81 and 2.00 depending on the numerical algorithm 
and the selection of collision partners (Meakin 1984). Changing the 
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sticking probability docs not dramatically change D„i- Collisions of 
polydisperse clusters always result in somewhat higher fractal dimen- 
sions than collisions between monodispersc clusters (Brown and Ball 
1985). Ballistic particle-cluster aggregation, BPCA, in which clusters 
grow by the addition of individual particles, give values of Dm close to 
3.0. The dimension D^ must always be smaller than Dm and is often 
assumed to be close to 2. This means that BCCA clusters have friction 
times that are nearly independent of radius. For such aggregates and 
typical conditions in the solar nebula at 1 AU, the relative velocities 
are between lO^'^ms^^ - Brownian motion - and 10^^ m s^^ - sedi- 
mentation. For turbulence-driven BPCA growth, the mean mass of the 
mass distribution scales with the evolutionary time, t, as i^, whereas 
the scaling for BCCA particles is exponential with time (Mizuno et al. 
1988). 

Numerical simulations of the growth process driven by Brownian 
motion under conditions appropriate for protoplanetary accretion disks 
(Kempf et al. 1998) show only a small dispersion of the friction times 
leading to small relative velocities and slow growth. The numerical 
calculations also demonstrated that the friction time increases rather 
slowly with radius: r oc i?"'^; Dm ^ 1.8, Do- ~ 1-6. No such simula- 
tions exist for turbulence-driven growth or coagulation driven by drift 
velocities. Experiments on the aggregation of micron-sized dust grains 
in a turbulent gas (collisional velocities in the experiment ~ 20cms^^) 
and during sedimentation (collisional velocities ~ 1 cms^^) confirm the 
formation of fractal dust grains with Dm between 1.9 and 1.7 in the two 
regimes (Blum et al. 1998 and references therein). The slow increase 
of Tf means grain growth is also slow. 

The friction time, Ty, must increase markedly {Dm > 2) to produce 
efficient, turbulence-driven coagulation (t/ > tg) (see equation (5)). 
Rapid growth will also occur when there is a reservoir of small particles 
(e.g. particles accreted from the molecular cloud core) leading to a 
BPCA-type growth process or during "runaway" growth where very 
large particles (decoupled from the mass distribution) sweep up a large 
amount of small grains (Safronov 1972; Weidenschilling and Cuzzi 1993; 
Schmitt et al. 1997; Weidenschilling 1997). In this case, the statistical 
coagulation equation cannot be used any longer (see, e.g., Tanaka and 
Nakazawa 1994). 

Experimental results suggest the possibility of enhanced coagula- 
tion rates of small iron particles by the presence of a magnetic field 
(Nuth et al. 1994; Nuth and Wilkinson 1995). Whenever a magnetic 
field of about 100 Gauss was present, the iron particles coagulated very 
rapidly and formed complicated aggregates. The relevance of this pro- 
cess for protoplanetary accretion disks depends on the possibility of the 
magnetization of small iron grains. 
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C. Critical velocities for grain sticking 

The growth process is governed by the sticking probability, S, and the 
strength of the adhesion forces holding the aggregates together. It is 
generally assumed that van der Waals forces are responsible for the 
sticking of small particles. The sticking probability for two colliding 
particles depends on the collision velocity, masses, shapes, and material 
properties. There may also be a critical velocity for sticking, Vcr- The 
existence of a critical velocity implies there is a steep transition from the 
velocity range where sticking occurs {v < Vcr, S ^ 1) to the velocity 
range where the particles bounce (v > Vcr, S ^ 0). Even for the 
simplest case of spherical particles, the calculated and measured values 
of Vcr disagree. There are many complicating factors to be understood 
in these collisions, so this disagreement is not surprising. 

Chokshi et al. (1993; see also Dominik and Tielens 1997) modelled 
the dynamics of sticking collisions between dust particles. Two colliding 
particles form a contact "neck", and the released energy due to the 
decrease of surface energy accelerates the particles towards each other. 
Elastic compression decelerates the particles, reverses the motion, and 
finally leads to acceleration in the opposite direction. When the grains 
again reach the point of first contact, they are still bound by a contact 
neck. After further separation, the binding neck becomes unstable and 
the contact can rupture. Whether or not rupture occurs depends on the 
amount of kinetic energy dissipated during the collision. In the model 
developed by Chokshi et al. (1993), energy dissipation is assumed to be 
dominated by elastic surface waves. The critical velocity they derived 
for central collisions between two spherical grains is: 

^5/6 
^- = l-0^ j;l/3^5/6^1/2 ' (8) 

where 7 is the surface energy per unit area, r is the reduced radius (r — 
'''7'^ , ri and r2 are the radii of the colliding particles), and p is the 
mass density of the spheres. i<^ is a material quantity which depends 
on the Poissons's ratios and Young's moduli and reflects the elastic 
properties. I A compilation of material properties for astrophysically 
relevant sytems is given by Dominik and Tielens (1997). For micron- 
sized quartz spheres, the critical velocity is calculated to be '^ 5 cms^^. 
Poppe et al. (1997, 1998) imaged the trajectories of micron- and 
sub- micron-sized particles colliding with a target to determine the stick- 
ing characteristics for different materials. For spherical silica parti- 
cles with smooth surfaces and with r = 0.6 /im, the critical velocity 
is 1.2 ms^^ which is an order of magnitude higher than predicted by 



f The factor in (8) has to be 1.07 instead of the 3.86 erroneously given in 
Chokshi et al., and the functional dependence has to be E^^'^ instead of E as 
misprinted in Dominik and Tielens. 
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equation (8). The shape of the particles is an important property for 
sticking. The sticking probabiUty increases significantly for irregularly 
shaped dust grains especially for submicron grains. In some experi- 
ments, no critical velocity is seen, and a significant fraction of irregular 
particles stick even at impact velocities exceeding 20 m s^^. 

Measurements with an atomic force microscope demonstrate that 
the pull-off forces are in agreement with the theoretical predictions. 
The discrepancy between experimental and theoretical values for Vcr 
in the case of the spherical particles is probably caused by an under- 
estimate of the actual energy dissipation rate in the theoretical calcu- 
lations. The general conclusion from the experiments is that micron- 
and sub-micron-sized grains stick with high probability as long as the 
relative velocities are below a few ms^^ - a limit which should not be 
reached by such particles in protoplanetary disks with a values below 
10-1. 

An interesting discovery of the experimental work is that coUisional 
charging occurs. The charge transfer in a coUison between micrometer- 
sized silicon particles and a flat quartz target surface results in averaged 
values of separated elementary charges per impact energy of 10^^ to 
10"^ C J^i. The charging is approximately proportional to the collision 
energy. Therefore, electrostatic forces between charged grains may be 
important for the growth process, but the sign and magnitude of the 
charges are difficult to predict a priori. 

Bridges et al. (1996) and Supulver et al. (1997) addressed the fur- 
ther growth from centimer-sized particles to meter-sized bodies. They 
found that several types of water frost-coated surfaces stick together 
when brought into contact at low enough sticking velocity, < 0.4 cm s^^ . 
Their experiments show that the structure of frost is important for 
sticking. They concluded that the formation of surface layers of frost is 
a necessary step for providing a sticking mechanism for the growth of 
larger particles in low-temperature regions of protoplanetary accretion 
disks. It might be interesting to note that the formation of frosty sur- 
faces plausibly results from sweepup of tiny grains onto the surfaces of 
the larger particles in relative motion. The layer would be regenerated 
following compacting bounces. 

If larger particles (cm- to meter-sized) reach relative velocities of 
the order 10 ms^^ or more, the collisons may cause compaction and 
disruption. High velocities for the larger particles would lead to frag- 
mentation and the production of a reservoir of small particles. 

D. Restructuring, compaction, and disruption in aggregate- 
aggregate collisions 

Aggregates with open BCCA-type structures cannot grow indefinitely 
through collisions. Collisions with larger grains or other aggregates 
should lead to restructuring and compaction at a certain mass and 
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velocity threshold. Compaction is necessary to increase the fractal di- 
mension (Dm > 2) which increases the frictional coupling and leads 
to a more efficient growth regime. This was the motivation for Wei- 
dcnschilling (see, e.g. Weidenschilling and Cuzzi 1993) to start with a 
fractal dimension of 2.11 in the size range 1 /im to 0.1 mm and then to 
introduce compaction starting with grain sizes of 1 cm. Although such 
a compaction is likely, the values of the size limits are unknown, and 
the change from Dm — 2.11 to 3 is arbitrary. And compaction may 
not occur only in collisions; recurrent heating events could also create 
compact particles. 

Theoretical calculations of the compaction by coUisons between ag- 
gregates by Dominik and Tielens (1997) predict a sequence of events 
depending on the ratio between impact energies and critical energies 
for the different processes: (i) sticking without restructuring, (ii) losing 
monomers, (iii) maximum compression, and (iv) catastrophic disrup- 
tion. Their evaluation of the different processes showed that rolling is 
the most efficient mechanism for restructuring. 

For velocities below 1 cms~^, the sticking probabilities in collisions 
of spherical Si02 aggregates are close to unity (Wurm and Blum 1998). 
There is no restructuring in turbulence-driven growth until coUisional 
velocities of 20cms^^. Measurements of the rolling friction between 
Si02 aggregates and a massive target by Wurm (1997) give a lower 
limit of 5.0 x 10~^''N, about a factor of 5 larger than the value used by 
Dominik and Tielens (1997). Bouncing and fragmentation dominated 
the results of collisions down to ^ Ims^^. Between 0.1 and Ims^^ 
(and especially close to the critical velocity) restructuring occured. 

The sequence found in the experiments is in qualitative agreement 
with the theoretical predictions of Dominik and Tielens (1997), but the 
threshold velocities for the different processes are higher by at least a 
factor of 2 than the calculated values, mainly because the calculations 
underestimate the frictional forces. These results imply that at least 
in the beginning of the sedimentation-driven growth, no restructuring 
will occur and the aggregates have to reach centimeter sizes before 
compaction occurs. 

IV. Observational consequences of particle size 

A. Theoretical changes of particle opacity 

As dust particles grow by collisions, the particle mass opacity, Hi, 
(cm^g^^), changes. Changes in opacity can in principle be observed 
in circumstellar disks, thus providing the signature needed to demon- 
strate that disks breed large particles and planets. The opacity at any 
wavelength. A, will not change until a significant fraction of the particle 
mass is in particles with sizes of order A or greater. 
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The flux density, F^,, from a disk viewed face on at distance, D, is: 

^- = hj ° B-^T{r)){l - e--'^W)27rrdr, (9) 

where the interior and exterior disk radii are r^ and i?£), B^ is the 
Planck function, T is the temperature, r^ is the optical depth, and r is 
the radial variable (cf Beckwith et al. 1990). The optical depth can be 
written in terms of the surface density, S, and mass opacity coefficient, 
K^ as: Ti,{r) = K^S(r). If the particles are small and the wavelengths 
are much larger than the size of the largest particles, the disk should be 
transparent, Ti, < \, and the emission is in the Rayleigh- Jeans regime, 
By « 2kTv'^. In this limit: 

F, ««,^2 / T{r)i:{r)rdr. (10) 

Regardless of the disk structure, the flux density is directly proportional 
to the particle mass opacity, k^. Although the absolute value of k^ is 
difficult to determine without a knowledge of T{r), S(r), and Rn [tq 
is usually not relevant in this limit), the frequency dependence can be 
observed directly. 

Particle emission in most disks becomes optically thin at millimeter 
wavelengths. The emission emerges from the outer parts of disks where 
T is a few tens of Kelvin, and S should be small, so that the approxima- 
tions made in Eqn. 10 are valid (Beckwith et al. 1990). Observations of 
disks over a range of wavelengths, between 0.6 and 2.7 mm, say, should 
provide a direct measure of the emissivity exponent, [3. If a mixture 
of particle types and sizes is present, the observed value is an average 
over the different constitutents. 

At long wavelengths, A > 0.1mm, k^ is expected to scale as a 
power of the frequency: k^ oc v^ . For compact spherical particles 
smaller than the observing wavelength, (3 — 2 for metals and insulators 
under a wide range of conditions (Bohren and Huffman 1983; Emerson 
1988). The power law exponent, (3, can be as small as 1 for certain 
types of materials - amorphous carbonaceous material, for example 
- and even smaller over limited ranges, but it must go to 2 at long 
enough wavelengths to fulfill causality relations, and it is expected that 
/3 « 2 at millimeter wavelengths for interstellar dust particles. Careful 
attempts to measure both k^ and (3 in interstellar clouds have, indeed, 
yielded /3 « 2 and values of k,, = (0.002--0.004)(A/1.3mm)-2 cm^ g-i 
(Hildebrand 1983; Draine and Lee 1984). 

Rocks, asteroids and planets are opaque to radiation at A ^^ 1 mm, 
in which case j3 = Q. Pebbles with sizes of order 1 mm should have 
exponents with intermediate values: < /3 < 2. If dust grains grow 
large enough to put most of the mass in bodies larger than pebbles, 
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it should have an observable signature at these wavelengths. Because 
the radiation docs not penetrate far below the surface of a large body, 
meaning most of the material never interacts with it, the absolute value 
of the opacity per unit mass must decrease as the particles grow. 

Calculations of (3 and Ki, for spherical particles show that they 
decrease as the particles grow (Miyake and Nakagawa 1993; Kriigel & 
Siebenmorgen 1994; Henning et al. 1995). Figure 3 illustrates the mass 
opacity coefficients, K^(a), for spheres of radii distributed as a power 
law with different maximum sizes, a^ax- In addition to the resonant 
absorption by silicates at 10 and 20/im, ice at 3/im, and possibly other 
solid-state features, the opacity is a power-law with l3 — 2 for A > 
100 /im. For X» a (Rayleigh scattering), Ki,{a) is independent of size 
a, whereas for A << a (i.e., in the geometrical cross-section regime), 
Ki/(a) is independent of v but has 1/a-dependence. For A 2± a, K^{a) 
is enhanced in comparison with an extrapolation from the Rayleigh- 
scattering regime, and Ki,{a) deviates from a power-law. Including the 
contribution of large particles with a ~ 1 mm or larger, Miyake and 
Nakagawa (1993) shows that /3 < 1 at millimeter wavelengths. The 
absolute value of Ki, can also be enhanced by particles with sizes similar 
to the observing wavelength. 
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Figure 3. The left hand side shows how the mass opacity coefficient, k^, 
varies as a function of the maximum radius, Umax , for a collection of spherical 
particles, whose number density, n, is a power law in particle radius, a, 
with a maximum radius, ttmax- The changing slopes at long wavelengths 
is equivalent to a change in /3. The right hand side shows the absolute 
value of Ki, at a wavelength of 1 mm as the maximum particle radius changes 
from dust (0.01 /xm) to boulders (100 m). After Figures 5 & 6 of Miyake & 
Nakagawa (1993). 



Pollack et al. (1985, 1994) proposed a model of dust particles which 
includes troilite (FeS) and some organics as major species in addition 
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to silicates, water-ice, and metallic iron. They found variations in the 
behaviour of k^, in the millimeter range, depending on the contributions 
of different constituents; the opacity is generally dominated by silicates 
but troilite and organic components can be more important than the 
silicates over restricted wavelength ranges. Opacity enhancements can 
also occur in specific materials as the result of lattice resonances. Since 
both changes in the particle chemistry and particle size can produce 
observable changes in the bulk opacities, it is not easy to demonstrate 
particle growth simply on the basis of changes in the opacities (Os- 
senkopf & Henning 1994, Henning et al. 1995). 



B. Laboratory measurements of particle opacity 

The interpretation of observations of circumstellar disks depends on 
a knowledge of the mass absorption coefficients of particle materials 
at the appropriate wavelengths, generally sub-millimeter to centime- 
ter. It is only recently that laboratory groups have measured optical 
data for cosmic dust analogues at the temperatures and wavelengths 
appropriate for circumstellar disks. Agladze et al. (1995, 1996) mea- 
sured absorption spectra of crystalline enstatite and forsterite grains 
and amorphous silicate grains with a ~ 0.1 — 1 fiia at T = 1.2 — 30 K for 
A = 0.7 — 2.9 mm and found that the mass opacity coefficients, K,y, for 
amorphous silicates were up to about ten times the values that are usu- 
ally adopted for interstellar dust particles (Draine and Lee 1984). The 
power-law behaviour also varies with T: /3 varies between 1.5 - 2.5 for 
the amorphous 2MgO-Si02 and MgO-Si02 but has a nearly constant 
value of 1.2 independent of temperature for amorphous MgO-2Si02. 
This behaviour comes from resonant absorption by low-lying two-level 
lattice states whose excitation depends on the temperature; such states 
are overlooked in the classic theory of solids that predicts j3 = 2. Thus, 
in addition to the other complicating factors determining mass opac- 
ity - size, shape, clustering effects, the presence of mantles, etc. — 
the material structure is enough to alter the observable properties. 
This result will complicate the interpretation of observations. Data 
on amorphous and crystalline silicates (Jager et al. 1994, Dorschner et 
al. 1995, Mutschke et al. 1998, Jager et al. 1998), oxides (Henning 
et al. 1995), sulfides (Begemann et al. 1994), and ices (Hudgins et 
al. 1993, Preibisch et al. 1993, Ehrenfreund et al. 1997) are of spe- 
cial importance for protoplanetary accretion disks (see Henning and 
Mutschke 1997 and Mennella et al. 1998 for a discussion of the tem- 
perature effects). Electronic data bases exist for ices and refractory 
solids and can be accessed by 



http://www.astro.uni-jena.de (refractory solids) 



http://www.strw. leidenuniv.n] (ices) and 
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V. Do we see particle growth in disks? 

Circumstcllar disks should contain large particles and planetesimals. 
There is little doubt that many observed disks are a few Myr old, and 
it should take only a few percent of this time to grow rocks from tiny 
dust grains. By a few Myr, much of the solid matter could be tied up 
in rocks and larger bodies, and we have seen that a shift in the mean 
particle sizes will have observable consequences for the spectral energy 
distributions. Do we see these effects? 

Changes in the spectral energy distributions (SEDs) of disks rela- 
tive to those typical of interstellar clouds are observed (Beckwith and 
Sargent 1991; Mannings and Emerson 1994; Koerncr et al. 1995). The 
flux densities, F^, of disks tend to fall more slowly at wavelengths 
longer than about 400 /xm than those of interstellar clouds, with the 
spectral index, a - where Fi, oa ly" - typically between 2 and 3, whereas 
a ~ 4 for the interstellar medium. Assuming optically thin emission 
in the Rayleigh- Jeans limit, the exponent, /3, is often between and 1 
for disks compared to the theoretical value of 2 for small particles in 
the long wavelength limit. The right hand part of Figure 4 shows the 
distribution of emissivity exponents, /3, determined from the observed 
spectral slope of disk emission near A 1 mm. Many of the disks in this 
sample have /3 < 1. It is tempting to interpret this change as the result 
of particle growth. 

Unfortunately, there are many adjustable parameters in the disk 
models, and several different parameter combinations usually produce 
acceptable fits to the same SEDs (Thamni, Steinacker, and Henning 
1994). A model typically has the following parameters: tq, Rjj, Ki,, 9 
(inclination to line of sight), Tg and q {T{r) = To(r/ro)^''), Sg and p 
(I](r) = So(r/ro)~^). There are rarely more than ten data points in a 
SED to be fitted. The SEDs are degenerate with respect to the models 
and almost certainly inadequate by themselves to demonstrate grain 
growth in disks. 

Fortunately, some of these parameters affect the millimeter-wave- 
length emission weakly or can be determined reasonably well from other 
wavelengths. The inner radius of the disk, tq, is unimportant for wave- 
lengths beyond 20 /im in almost all cases. Both Tq and q are well 
constrained by the optically thick emission between 10 and 100 /im. 
The disk inclination, 9, is unimportant for the optically thin (mm-A) 
emission but important for the optically thick (A < 100 /im) part of the 
SED. In some cases, 6 can be estimated directly from images. 

The remaining parameters - Eg, p, Rd, Hv ~ play a role in the 
long wavelength emission and, hence, the derivation of the emissivity 
exponent, /3, from the observed spectral slope, a. In the simplest case of 
p = (constant surface density), t^ < 1 for all wavelengths longer than 
100 /im, say, and T{r) large enough for all r to make the Rayleigh-Jeans 
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approximation valid for millimeter wavelengths, /3 = a — 2 (see Eqn. 
10), the spectral index gives the emissivity exponent directly. But if the 
density decreases at least as fast as r~2 (p = |)^ there is substantial 
contribution to the long wavelength emission from the inner parts of 
the disk, li p > |, the inner parts of the disk dominate the emission 
making R^ irrelevant but implying that most of the emission at any 
wavelength comes from radii near the radius at which the disk becomes 
optically thick (see Beckwith et al. 1990 for a detailed discussion of 
these effects). In this case, the spectral index, a, depends on (3 and 
both p and q. Specifically, Beckwith et al. (1990) derived /3 in terms 
of the model parameters as: 

/3«(a-2)(l + ^r ^ J , (11) 



(2 - q) In 



(2-p)t 



where t,^ — ^^"^^"^ is the average optical depth at frequency, i^, and 
Mu is the mass of the disk. J Equation 11 depends strongly on p which 
cannot be determined from SEDs alone. For this reason, the emissivity 
exponent derived from and SED is model dependent. 

The right hand part of Figure 4 includes the distribution of /?'s 
derived from model fits assuming the surface density is proportional 
to r^^/^. Including the surface density distribution tends to increase 
the derived value of the emissivity exponent relative to the observed 
spectral slope. For example for 1.3 mm observation of a disk with Mo = 
0.03 Mq, Rd = lOOAU, and k^ = 0.02cm2g"\ the average optical 
depth, f^ — 0.24, and /3 w 1.36(q; — 2). Any conclusions about particle 
growth depend on knowledge of how much optically thick parts of the 
disks contribute to the spectral energy distribution at long wavelengths. 
It is probably for this reason that the samples of disks with measured 
a's at long wavelengths have not be enlarged much since the work of 
Beckwith & Sargent (1991) and Mannings & Emerson (1994). 

To measure the optical depth, it is essential that the disk emission 
be resolved at wavelengths where it is relatively transparent so that 
the surface density distribution may be determined directly from the 
distribution of optical depth. Resolving the disks is time-consuming 
with the present suite of millimeter facilities. Several groups have now 
managed to resolve the emission from one star, HL Tau, by carrying 
out novel experiments with existing millimeter-wave telescopes. Lay 
et al. (1994) and Lay, Carlstrom & Hills (1997) combined the CSO 
and JCMT into a single-baseline interferometer to observe at 0.65 and 
0.87mm and combined it with data from the Owens Valley millimeter 
interferometer at 1.3 mm. Mundy et al. (1996) extended the BIMA 
array to measure a size of the disk at 2.7 mm wavelength, and Wilner, 



J Beckwith et al. derive the spectral index from i^Fi,, not Fi, used here. 
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Ho & Rodriguez (1996) used the VLA to get sub-arcsecond resolution ol 
the disk at 7 mm. Each group resolved the emission but with beam sizes 
of the same order as the disk size, and the sub-millimeter observations 
had u — V plane coverage too poor to map the brightness distribution. 
All groups measured the average brightness temperature and overall 
disk orientation and constrained subsequent model fitting much more 
than is possible from the SEDs alone. 
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Figure 4. The left hand part shows the spectral energy distribution of HL 
Tau from about 350 /im to 6 cm wavelength. The solid line fits the short 
wavelengths - thermal emission by dust - and and the dashed line fits the 
long wavelengths - free-free emission from ionized gas. The long wavelength 
spectral index of the dust is 2.9, meaning j3 — 0.9 with no corrections for 
optically thick emission (Wilner, Ho, & Rodriguez 1996). The right hand 
portion shows the distribution of /3's derived in in two ways for a larger 
sample of disks. The "observed" /3's are derived directly from the spectral 
index as described in the text: Pots —a — 2. The "model" /3's are derived by 
assuming a surface density distribution as discussed in the text (Beckwith & 
Sargent 1991). 

Unfortunately, the combined data are not consistent with simple 
power-law models for disk optical depth. The apparent size of the 
HL Tau disk is approximately constant from 0.65 mm to 7 mm. Be- 
cause the optical depths should drop by factors of 10 (for (3 — 1) to 
100 {/3 — 2) over this wavelength range, any centrally concentrated 
density distribution, p > I, would immediately imply smaller sizes at 
longer wavelengths. Lay, Carlstrom & Hills (1997) discuss the modi- 
fications in the models that might reconcile the data. They conclude 
that a mixture of particle sizes such that there is no unique value of 
P may explain the data. We note that there would have to be sub- 
stantial contributions to the emission from particles of order 1 cm for 
this conclusion to hold. Although such particles are not yet the size 
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of planetesimals, they are more than four orders of magnitude larger 
than typical interstellar grains and would be present only if there was 
substantial particle growth in the disk. HL Tau is one of the youngest 
disks (Cohen 1983; Beckwith et al. 1990), so demonstration of parti- 
cle growth in this system would imply that it occurs quickly as theory 
suggests. 

HL Tau is also one of the more complicated star/disk systems to 
interpret, owing to the large amount of surrounding material and the 
likelihood that high optical depths arc contributing to the millimeter- 
wave emission (cf. Stapelfeldt et al. 1995; Beckwith and Birk 1995; 
Weintraub et al. 1995; Cabrit et al. 1996; Menshchikov et al. 1998). 
Mundy et al. (1996) observe a brightness temperature of 29 K from 
the HL Tau disk at 2.7 mm with an outer radius of order 150 AU. This 
temperature is very close to the physical temperature of the outer re- 
gions of the disk: Beckwith et al. (1990) derive T = 32K at a disk 
radius of 100 AU. Since the brightness temperature, Tb, is related to 
the physical temperature by Tb = T{1 — e""^), there must be a sub- 
stantial contribution from optically thick dust that is not easily taken 
into account in the model fits. 

Because HL Tau is much younger, more massive, and considerably 
more complex than most of the disks in Figure 4, the derived /3's may 
still be valid in many cases. But since the first direct measurement of 
the optical depth distribution show the model assumptions are invalid, 
we cannot conclude yet that the long wavelength spectral indices tell 
us much about the particle properties. 

The evidence for grain growth is, therefore, inconclusive. There 
are quite a few disks with spectral indices that are inconsistent with 
optically thin emission from very small grains. The spectral indices are 
not unambiguous indicators of large grains, however. Contributions 
from high optical depth emission within the disks can flatten the SEDs 
in similar ways to changes in the grain-size distribution. Furthermore, 
/3 = 1 is not sufhcient to demonstrate grain growth. Certain kinds of 
materials can produce an index of 1 even for small grains. There are 
quite a few disks for which the spectral index implies /3 < 1, but none 
have been resolved at millimeter-wavelengths. 

These observations are difficult, because the required spatial reso- 
lution is much higher than that routinely available with existing milli- 
meter-wave interferometers. This situation will change during the next 
decade. The proposed construction of a large millimeter-wave interfer- 
ometer in the southern hemisphere - the American MMA and European 
LSA projects - will yield a combination of sensitivity, resolution, and 
wavelength coverage allowing all the disks in the dark clouds within 
< 300 pc to be easily resolved. Mapping the disk emission should 
provide a direct measure of the optical depth, assuming the physical 
temperature can be observed independently, from molecular lines, say. 
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From the optical depth at different wavelengths, it should be possible to 
demonstrate that a population of large grains or rocks does contribute 
to the emission, if such a population exists. This avenue is the most 
obvious hope for a demonstration of grain growth in disks. 

In the context of planet formation, it is gratifying that the evidence 
favors grain growth in some disks, even if this evidence is insufficient 
to establish unambiguously the presence of large particles. The situ- 
ation is similar to that of the disks themselves in the last decade. At 
that time, many observations were consistent with the existence of cir- 
cumstellar disks with the right properties for planet formation, but the 
evidence was not insufficient for proof beyond reasonable doubt. Now, 
there is no doubt that disks exist with the general properties then in- 
ferred, and a separate line of research has uncovered planets themselves 
- albeit indirectly - orbiting other stars. It is natural to suppose that 
disks are the necessary precursors to planets. The next generation of 
millimeter-wave observatories may well prove the evolutionary link by 
demonstrating particle growth, too. 
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